On the maximum drawdown during speculative 

bubbles 

Giulia Rotundo 1 ' 2 * Mauro Navarra 2 
1 Faculty of Economics, University of Tuscia, Viterbo, Italy 
2 Faculty of Economics, University of Rome "La Sapienza" 

February 2, 2008 



Abstract 

A taxonomy of large financial crashes proposed in the literature locates 
the burst of speculative bubbles due to endogenous causes in the frame- 
work of extreme stock market crashes, defined as falls of market prices 
that are outlier with respect to the bulk of drawdown price movement 
distribution. This paper goes on deeper in the analysis providing a fur- 
ther characterization of the rising part of such selected bubbles through 
the examination of drawdown and maximum drawdown movement of in- 
dices prices. The analysis of drawdown duration is also performed and it 
is the core of the risk measure estimated here. 

Key words Risk measure, drawdown, speculative bubbles. 
PACS 89.65.Gh Economics, business, and financial markets, 
89.90.+n Other areas of general interest to physicist 



1 Introduction 

An insight in the long term behavior of portfolios is a delicate task in long term 
investment strategies. The need to consider extreme financial market events 
encompasses the investigation of large financial crashes, that were already clas- 
sified as outliers with respect to the bulk of market drops, and often associated 
to the burst of speculative bubbles due to endogenous causes [T]. 
This paper aims at extracting risk features that characterize the rising part of 
such speculative bubbles. 

Our data selection relies on the huge data analyses worked out by Johansen 
and Sornette [U H] ■ I n their several papers, they develop and support through 
empirical evidence the theory describing speculative bubbles due to endogenous 
causes like as systems close to some rupture point. In particular large market 
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indices drops ending speculative bubbles due to endogenous causes have been 
characterized through the occurrence of log-periodic-power laws (LPPL), in- 
terpreted as a signature of underlying cooperative phenomena among market 
agents (Fig. [TJ Tab. [JJ). These results were related to analyses of drawdown 
movements of market indices prices (DD), thus providing systematic taxonomy 
of crashes [TJ [3J H] ■ Our analysis starts from this point and goes on inquiring 
drawdown market prices movements duration and size. The analysis of dura- 
tion of drawdown movement is relevant in itself, considering the key role of time 
in catastrophic events like as large financial crashes are, that show a drop of 
prices within short time intervals [TJ (3J [5] . The analysis of drawdown size re- 
quires as a preliminary step the selection among the several different definitions 
of drawdown that were proposed in literature. Results depend strictly on the 
measure chosen for drawdown, and on the actual data set [TJ [5j [6J [7] . A first 
analysis classified large stock markets indices crashes as outliers whether con- 
sidering the stretched exponential distribution of market drops 6 . A further 
analysis for noise filtered drawdowns allowed to calibrate the parameter in order 
to avoid outliers [5] . Other studies using different definitions of drawdowns and 
considering different data sets proposed for returns drawdown probability den- 
sity fitting stable models [5], the generalized Pareto distribution (GPD) pJ[H], 
and the modified GPD. The worst scenario of maximum drawdown movements 
(MDD) has also been examined. Results depend also on the time window on 
which the maximum is estimated. For example in [7] the time window selection 
corresponds approximately to half an year, the entire data sets cover more than 
20 years, and the results are used for probability distribution calibration. In 
[8] moving windows of 63 data each are used on datasets including at least 15 
years. In both cases the long time series of market indices data that were used 
include lower and bigger financial crashes, the rise and the deflate of bubbles, 
till to the return to fundamentals, and the several changes of regimes included 
in societal and trading dynamics spanning the last 15-20 years. Therefore, the 
maximum drawdown measure calculated on moving windows gives samples of 
a process obeying several different dynamics. Here we aim at using a different 
approach for drawdown modeling. Instead of splitting the entire time series that 
are available, we extrapolate the probability of drawdown and maximum draw- 
down during the rising part of the speculative bubble, looking for the extraction 
of common features. 

We compare the results obtained on the DD through the estimate of risk mea- 
sures. Risk exposure in financial markets has been described through several 
measures (VaR, CVaR, etc.) based on statistical properties of data that do not 
consider the order of data sequences. In fact empirical mean, variance, as well 
as higher order moments of probability distributions are invariant under data 
shuffle. But long downward trends containing long lasting sequences of consec- 
utive drawdown price movements could suggest investors to withdraw from the 
market, and they can quite force small investors to such a choice [5]. There- 
fore, measures of risk based both on the duration of consecutive market drops 
and on the maximum drawdown can play a relevant role in driving investment 
strategies. 
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We stress again that we use here a methodological approach deeply different 
form other comparative studies on risk measures, because we aim at extracting 
features of periods (the rising part of speculative bubbles due to endogenous 
causes), that are supposed to be driven by the same kind of dynamics, and we 
proceed comparing them with the entire data set available, whilst other studies 
that compare risk measures on long time windows do not take into account the 
evolution of dynamics driving stock markets across decades [8] . 

2 Large stock market indices crashes 

A series of papers about speculative bubbles proposes an explanation of large 
financial crashes due to causes endogenous to the market. Similarities with 
critical phenomena like earthquakes and the sound emission in materials close 
to the rupture point led the research to the detection of cooperative underlying 
phenomena evidenced through discrete scale invariancc in financial data (FX, 
Gold, stock market indices) close to large crashes that can thus be well described 
as critical points 2J . In the case of market indices the logarithm of index values 
close to the crash time is described by a characteristic LPPL: 

F(t) ~ A + B{t c - t) m + C{t c - t) m cos{ujln(t c -t)-<j>) t < t c (1) 

where t c is the crash time, and A, B, C, m , to, </>, and t c are parameters to be 
estimated via numerical optimization. A further taxonomy that distinguishes 
the causes that generate crashes can be found in pQ and confirms the detection 
of LPPL as a hallmark of speculative bubbles due to causes endogenous to the 
market [U QUI [11] cross-validated with the amplitude of the crash ending the 
bubble as an outlier from probability distribution of the bulk of price variations 
[4j [6] . Stock market indices datasets reported in Tab. Q] were selected lying 
on the results of pQ in which the selection of an outlier market drop is used in 
itself as a definition of crash and then it is cross-validated with the occurrence 
of LPPL. The rising period of the bubble in which the LPPLs occur is selected 
starting from a lowest point before the rise of the bubble and ending at the 
maximum value of the index [5] . 

3 Drawdown duration distribution 

The time plays a key role in crash definition [TJE] and in the reaction of investors 
to the change of market condition, but most attention up to now was paid to 
the MDD than to the time in which it occurs, although the time measure was 
proposed as a first step towards the definition of some "drawdown velocity" 
[1 . Therefore, for each selected data set we exploit the duration of drawdowns, 
defined in the most simple way as the number of time steps from a local peak 
to the next local minimum, corresponding to the so called "pure" drawdowns 

of m Ei Ei . 
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A first empirical result for drawdown duration is shown in [7] , where the neg- 
ative binomial or the Poisson distribution is proposed for modeling drawdown 
duration, and the Gamma or Pareto distribution for modeling the maximum 
drawdown. 

Let {pt}t=i.T be a time series of a stock market index daily closure price, and 
let r t = log(p t /pt-i) be the usual log-returns. The following definitions are in 
accord with the drawdown definition used in [7J, apart from a normalization 
factor that is not relevant for our analysis (Fig. [2]). 

Definition 1 Let pk be a local maximum, and pi be the next local minimum: 
p k > Pk+i > ■ ■■ > pi, with I - k > 1, k > l,andp fc > p k -i, Pi+i > Pi- (2) 
A pure DD (size) is defined as the partial sum of daily returns 

l-k 

DD k .i =\ r k +i + r k +2 H n \=\ ^ r k+] \=\ log— | (3) 

Definition 2 The drawdown duration is the time length I — k of the sequence 
of negative returns defined through . 

DD duration was estimated on each selected time series. The power law was 
fitted to DD duration histograms. On the sets of the rising parts of selected 
bubbles shown in Fig. [T]the power law exponent ranges from 0.46 to 1.7, having 
a mean value a ~ 1. On the entire data sets a ranges from 0.67 to 1.10, with a 
mean value a ~ 0.94, accordingly to the presence of longer and deeper crashes 
in the enlarged time domain. The Poisson distribution gives a poor fit of DD 
duration histogram on the lowest duration values, and no useful result comes 
from the negative binomial fit. Therefore, we are going to use the empirical 
estimates for the risk measure computation (Fig. [3]). 

Studies on speculative bubbles have attempted to extract common features 
beyond LPPL. A first attempt [7J to relate GARCH parameters with the pre- 
diction of the volatility connected to the drawdown movement is made clearer 
in [12]: the presence of a bubble drives the high volatility in the GARCH model, 
but the reversal implication does not hold. The power law exponents of draw- 
down duration decay seem to evidence some general feature of drawdown dura- 
tion interesting in itself. 

The grouping of the exponents in a range supports the hypothesis of close risk 
profiles across international markets, form the DD duration point of view. How- 
ever, a further correlation analysis comparing the a estimated on each series and 
the LPPL paramenters oj and z does not show significative empirical correlation, 
like at it was stated for GARCH parameters. 

4 Relationships with other drawdown definitions 

The delicate task of defining what is a crash, a large drop or a significant change 
of regime is far from being well assessed [12] . Also the closely connected defini- 
tions of drawdown market movements analyzed in the literature are not homo- 
geneous [TJ [9l [13] . Therefore we are going to examine the relationship between 
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the so called pure drawdown (DD), e-drawdown (eDD) [HE], that are used in 
the theory of speculative bubbles, and the maximum drawdown (D) used in 
[HI H3]j that provides a risk measure interestingly related to the Calmar and 
Sharpe ratio. 

Definition 3 The maximum drawdown M is the longest partial sum of daily 
returns 

M = max{DD k)l | conditions © hold} (4) 

The definition of eDD relaxes the condition @ allowing for small rises of values 
during the decreasing sequence. This description can be formalized through the 
following set of conditions 

k = arg(max s (s - I)), s.t. p k > p k+1 > p k+2 - e > ■ ■ ■ > pi - e, I - k > 1, k > 1, 
p k > p k -i,p k >Pi,i= 1, ...,l,pi +l > pi + e,e > 0. 

(5) 

and therefore to the following 

Definition 4 The eDD ending at time I is given by 

eDDi = {DDi | conditions © hold}. (6) 

This definition states that price growth below a certain magnitude e is ignored, 
and it serves to filter noise. If e = the definition corresponds to "pure" 
drawdown; filter size was chosen as a function of the empirical volatility a 
(e = 0, er/4, a/2) [lj. This definition includes in the particular case e = the 
definition of pure drawdown ending at time t (DDt). 
Definition 5 The duration of eDDi is given by 

max s {(s ~ I) | conditions ([5]) hold} (7) 

D will be the worst loss, i.e. max of the above. 

The definition used in [13] starts from a stochastic processes approach: 
Definition 6 Let {X(t)} a stochastic process. D is given by 

D = sup te[0 ^ T] (sup se[0 ^X(s) - X(t)). (8) 

This measure gives the range from the maximum to the minimum anytime the 
maximum precedes the minimum. D was calculated using X(t) — log(p t ), so to 
deal with returns. 

Although the definition of eDD was aimed only at filtering the noise the following 
remark holds: 

Remark 1 Let X(s) = log(p s ). For each t ending time of an eDD t : 

eDD t < D 

Remark 2 D is an eDD t for X(t) = log(pt), e > maxt \ pt — Pt-i I- 
Henceforth risk measures based on DD and D will provide bounds to eDD. 
The following chain of inequalities is a direct consequence of the weakening of 
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conditions in the definitions. 

Lemma 1 In accord with the definitions reported above: 

DD t < eDD t < D 

M < D W 

The size of drawdowns described by D is higher of M, due to the weakening 
of the descent conditions, and provides a worse scenario. 

The same remark allow to conclude that the same inequality holds on the du- 
ration of DDt and eDDf . 



5 A measure of risk based on drawdown move- 
ments size and duration 

Stock market indices actually are a particular portfolio, basing on a weighted 
mean of selected stock prices, and to buy/sell stock market indices has the 
meaning to buy/sell a previous selected financial product replica of the index 
(Exchange Traded Funds, certificates). Portfolio risk measures considering DD 
and MDD should be used in a complementary way with respect to the tradi- 
tional ones (VaR, ES) at least in the case of Stable Paretian distribution b 8 . This 
section aims at comparing the behavior of a risk measure based on drawdown 
movement size on the rising part of speculative bubbles and on the entire time 
series. Results of [1] that classify as outliers with respect to the stretched expo- 
nential distribution the large financial crahes rising from the burst of specula- 
tive bubbles due to endogenous causes were extended considering coarse-grained 
drawdowns (e-drawdowns) [5] , and later the GPD distribution for negative tails 
was tested [7J |H1 HH US] . These last results go beyond the mere distribution hy- 
pothesis testing, relying on literature on extreme events, that proposes the GPD, 
as a universal description of the tail of distributions of Peaks-Over-Thresholds. 
The same approach can be used for the MDD estimates. Although the duration 
plays a key role in drawdowns little or no analysis has been reported on the 
literature on drawdown duration [7j. Therefore, we perform analyses on the 
joint probability of both drawdown size and drawdown duration. The following 
measure of drawdowns 

oo 

Pr{DD < s} = Pr{DD k j <s \ l-k = d}Pr{l - k = d} (10) 

d=l 

is estimated on the rising part of each selected bubble and on the entire time 
series of the corresponding stock market index. This approach differs from the 
one used in [7J [H] , that relies on probability obtained as a best fit of empirical 
drawdown size distribution, but that is based on independence hypothesis for 
drawdown duration modeling. 

The weight of the biggest crashes in the entire time series evidences in lowering 
curves corresponding to entire time series. 
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Following the approach of [7] we fit the GPD 

G tf) (x) = l-(l+S*x/0)? , xeD(£,p) (11) 

on the set of (fTO]) estimated on the rising part and on the set of entire time 
series (Figs [4] [5]). This result provides the mean behavior of drawdown to 
which refer in case of speculative bubbles, and evidences the difference with 
respect to the common behavior of the entire time series. On the rising part 
of the bubble (3 = 0.0089(0.0085,0.0093), £ = 0.4273(0.3783,0.4763); on the 
entire time series /3 = 0.0153(0.0148, 0.0158), £ = 0.3012(0.2698, 0.3326). Values 
reported in parenthesis are 95% confidence interval. The ranges of £ and [3 
parameters are not overlapping. The change of parameters value gives also 
another measure of large drops considered in the entire time series if compared 
with the drops experienced during the bubble. For each selected bubble Fig. \5\ 
shows the cumulative function (flTJ|) of both the rising part of the bubble and the 
entire index. The closest the functions, the closest the drawdown structures are, 
and the lowest the impact of the outlier, that implies also that minor crashes 
in the rising part of the bubble have a relevant weight. As an example, drops 
during the rising part of the speculative bubble of Nasdaq collapsed in April 
2000 had an intensity close to the Nasdaq crashes of the previous years, leading 
to highly overlapping curves. 

Referring to Fig. [4] "fit 2" is the regression on the sets of (fT0|) estimated on the 
rising part of speculative bubbles, "fit 2" separates DD behavior on the entire 
data set from 1/3 of (JTTJJ) estimated on the rising part of speculative bubbles, 
corresponding to series exhibiting smaller and/or shorter DD. 
Risk profiles of rising part/entire time series are not completely separated, but 
the result allow to use lower risk profiles in the rising part of speculative bubbles 
whether carefully encapsulated in a predictive scheme of burst of bubbles like 

6 Maximum of drawdown movements 

Other measures of drawdown behavior were considered, in the case of the oc- 
currence of the worst scenario. A discussion of coherence properties of measures 
DD-based is carried on [3 [9]. Pr{M < s} was calculated in [7] conditioning 
to the duration, like in (jXOJ) . The probability distribution of M was extracted 
through a segmentation of 15-20 years long time series. This procedure relies 
on the (undeclared) assumption that data in time windows considered obey a 
homogeneous underlying process, that of course is a very raw assumption. The 
need to split windows is due to the fact that only one M is available for each 
time series, but several samplings are needed in order to build a distribution. 
Here we follow a deeply different methodological approach, aiming at extracting 
features that can be shared by the rising part of speculative bubbles. There- 
fore we consider max(DDk^) conditioned to the duration and we estimate and 
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compare 

oo 

E[M] = max{DD k j \ l-k = d}Pr{l - k = d} (12) 

d=l 

on each time series. Results are shown in Fig. [6l The leftmost picture reports 
the results on the rising part of the bubble, the rightmost on the entire time 
series. This estimate compares the size of maximum sequences of drops inside 
the bubble with the ones outside it, allowing to understand the impact of a 
huge financial crash with respect to the smaller crashes internal to the bubble. 
E[M] £ (0.01,0.08) (rising part) and E[M] € (0.05,0.13) (entire series). Ranges 
are partially overlapping, remarking that some crashes classified as huge in some 
indices had a magnitude not so impressive for other indices. 
Referring to Fig. [SJ the discriminant analysis locates E[M] S (0.01,0.05) for 
the rising part of speculative bubbles and £7[M] € (0.08,0.13) for the entire 
time series. Intermediate values (0.05, 0.08) report the situation for the crashes 
of Nasdaq 100 in 1998 and 2000, DAX 40 in 1998, Argentina in 1997, Hong 
Kong in 1994 and 1997 (rising parts of the bubbles; on Fig. Q] labeled, re- 
spectively, Nasdaq 100_1998, Nasdaq 100_2000, DAX 40_1998, Arg MervaL1997, 
Hang Seng_1994 and Hang Seng_1997), and for FTSE 100, Chile General (Igpa), 
Jakarta SE Composite (entire time series). 

Fig. [7| shows D € (0.04,0.17) for the rising part of speculative bubbles; D 6 
(0.08, 0.35) for the entire time series. 

The weakening of the descent conditions can also be addressed for the higher 
overlap in windows values. Anyway this measure is worth of being considered 
because it gives the very worst loss case in case of long time buy-and-hold 
strategies. 

7 Conclusions and further developments 

This paper focuses on risk measures based on drawdowns and has introduced 
the empirical estimate of duration in the loss function underlying a risk measure. 
The analysis of the duration of drawdown is relevant for fund managers, that 
can bear high volatility periods, but that risk to lose clients whether a large 
sequence of drops happens. Empirical results are used to extract characteristics 
common to the rising part of speculative bubbles and to evidence similarities 
and differences with respect to the entire time series. A comparison between 
risk measure DD-based shows the improvement in the discriminant analysis 
obtained by the introduction of DD duration. The approach used here is deeply 
far from the most common usage to split long time series or to simulate them 
basing on the raw hypothesis of a homogeneous process underlying long time 
series. Recently an alert system |12j was proposed in order to open the way 
to practical usage of the LPPL bubble theory for the forecast of the bubble 
ending time. We look forward for the embedding the present analysis into the 
predictive scheme of [H] for the calibration of investment strategies during the 
rising part of bubbles, before the expected crash time. 
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FTSE 1 00 J 987 FTSE 1 00 J 997 FTSE 1 00 J 998 
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Kuala Lumpur SEEmasJ 994 Hang Seng J 994 xio 4 Hang Seng J 997 




92.80 93.20 93.61 94.00 92.02 92.68 93.34 94.00 95.37 96.11 96.86 97.60 



Figure 1: Rising part of speculative bubbles selected. Large financial crashes 
were chosen in accord to [T|. The complete list is given in Tab. [1] Time is 
reported on x axis, in accord with the notation used by jT]. y axis reports the 
value of each index. 
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Table 1: List of crashes. Data sets of rising part of speculative bubbles are 
chosen since the rise of the bubble to the expected crash time [2] 



Index 


Crash year 


Rising part 




Entire dataset 


Nasdaq 100 


1987 


10/01/1985 


- 10/02/1987 


1 n /m /1 Q87 


1 9 /^i /9nni 


Nasdaq 100 


1998 


01/02/1995 


- 07/20/1998 


1 n /m /1 Q87 


1 9 /^i /9nni 


Nasdaq 100 


2000 


04/01/1997 


- 03/27/2000 


10/01 /1987 - 


1 9 /31 /9001 


FTSE 100 


1987 


07/01/1985 


- 07/16/1987 


07/01 /1985 - 


19/31 /9001 


FTSE 100 


1997 


04/03/1995 


- 10/03/1997 


07/01/1985 - 


12/31/2001 


FTSE 100 


1998 


03/14/1995 


- 07/20/1998 


07/01/1985 - 


12/31/2001 


DAX 40 


1998 


08/06/1996 


- 07/20/1998 


08/06/1996 - 


12/31/2001 


Brazil Bovespa 


1997 


03/06/1996 


- 07/08/1997 


03/06/1996 - 


12/29/2000 


Arg Merval 


1997 


03/15/1995 


- 08/20/1997 


03/15/1995 - 


12/29/2000 


Venezuela SE Gen 


1997 


09/12/1995 


- 09/23/1997 


09/12/1995 - 


12/29/2000 


Chile General (Igpa) 


1991 


10/10/1990 


- 10/08/1991 


10/10/1990 - 


12/31/1996 


Chile General (Igpa) 


1993 


05/10/1993 


- 02/04/1994 


10/10/1990 - 


12/31/1996 


Mexico Ipc 


1994 


06/15/1993 


- 02/08/1994 


06/15/1993 - 


12/29/2000 


Mexico Ipc 


1997 


10/18/1995 


- 08/06/1997 


06/15/1993 - 


12/29/2000 


Kuala Lumpur SE Emas 


1994 


08/06/1992 


- 01/04/1994 


08/06/1992 - 


21/31/1996 


Hang Seng 


1994 


09/12/1991 


- 01/04/1994 


09/12/1991 - 


12/29/2000 


Hang Seng 


1997 


01/02/1995 


- 08/07/1997 


09/12/1991 - 


12/29/2000 



ED = | log (|>. t,l 




Figure 2: Pure DD definition in accord to ^ 
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Nasdaq 100 J 987 



Nasdaq 100 J 998 



Nasdaq 100_2000 




2 4 6 
Chile General (lgpa)_1991 



2 4 6 8 10 
Chile General (lgpa)_1 993 



2 4 

Jakarta SE Composite_1994 





2 4 6 8 10 
Kuala Lumpur SEEmasJ 994 



4-0 o 




10 2 4 



Figure 3: Histogram of the duration - rising parts of the speculative bubbles. 
Usual histogram bars were substituted by diamonds because of graphic clarity. 
x-axes report the DD duration in term of days; y-axes report their frequencies. 
On each example the solid line is the power law regression curve; the dotted line 
is drawn using the bounds of the 95% confidence intervals. The decay exponents 
a G [0.46, 1.7] where 0.46 is the minimum value and 1.7 is the maximum value 
over all the reported samples; the mean value of a is approximately at 1. 
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Figure 4: Best fit of empirical estimates of (fT0|) through a;-axis reports the 
values of s; y-axis the cumulative probability. The curve labeled "fitl" is the 
regression curve obtained by fitting (|lip on the empirical estimates of (fTU)) on 
the set of entire time series. Parameters values are (3 = 0.0153(0.0148,0.0158), 
£ = 0.3012(0.2698,0.3326). The curve labeled "fit2" is the regression curve 
obtained by fitting (fTTjl on the empirical estimates of (fT0|) on the set of 
rising part of bubbles. Parameters values are (3 — 0.0089(0.0085,0.0093), 
£ = 0.4273(0.3783,0.4763). Values reported in parentheses are 95% confidence 
interval. 





Figure 5: Empirical estimates of (fT0|) . rc-axes reports the value of s; y-ax.es the 
cumulative probability. Each figure reports (|10|) estimated on the rising part of 
a selected bubble (solid lines) and compared with lfT0|) estimated on the entire 
time series (dashed lines). 
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Rising pari ol [lis speculative bubble Emise time series 




Figure 6: Histogram of maxima of pure drawdown (|12p . x-axes reports values 
of E[M]; y-ax.es reports the counting of them. On the set of rising part of 
speculative bubbles the range is (0.01,0.08); on the set of entire time series the 
range is (0.05,0.13). 



Rising part of the speculative bubble Entire data sei 




Figure 7: Histogram of D ([8]). x-axes reports values of D; y-axes reports the 
counting of them. On the set of rising part of speculative bubbles the range is 
(0.04,0.17) ; on the set of entire time series the range is (0.08,0.35). 
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